In the Claims: 



1-2. (Cancelled) 

3. (Previously presented) The method of Claim 4, wherein located at 
each collocation point tj is a value of u(tj), respectively, to be interpolated with 
polynomials. 

4. (Currently Amended) A method of simulating a circuit, the method 
comprising: 

defining a differential-algebraic equation of the circuit; 

defining a simulation time interval corresponding to the differential- 
algebraic equation; 

dividing the simulation time interval into time intervals, wherein the time 
intervals have corresponding polynomials for each time interval, wherein each 
polynomial is a portion of an approximation to a desired solution of the 
differential-algebraic equation; and 

applying a collocation method to discretize the differential-algebraic 
equation; 

wherein: 

the simulation time interval has collocation points, and wherein the 
interpolating polynomial has a degree of M; 

the approximation to the desired solution of the differential-algebraic 

^(0 = Z«*r A (<), 

equations is *=° wherein Im represents an M-point interpolation 

operator, u(t) represents a solution, — Uk represents Fourier coefficients, Tk(t) 
represents the interpolating polynomial, and M Is represents the highest degree 
of the interpolating polynomials. 
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5. (Currently Amended) A method of simulating a circuit, the method 
comprising: 

defining a differential-algebraic equation of the circuit; 

defining a simulation time interval corresponding to the differential- 
algebraic equation; 

dividing the simulation time interval into time intervals, wherein the time 
intervals have corresponding polynomials for each time interval, wherein each 
polynomial is a portion of an approximation to a desired solution of the 
differential-algebraic equation; and 

applying a collocation method to discretize the differential-algebraic 
equation; 

wherein: 

the simulation time interval has collocation points, and wherein the 
interpolating polynomial has a degree of M; 

the approximation to the desired solution of the differential-algebraic 

M 

eiquations is k=0 wherein Im represents an M-point interpolation 

operator, u(t) represents a solution, -Uk represents Fourier coefficients, Tk(t) 
represents the interpolating polynomial, and M is represents the highest degree 
of the interpolating polynomials; and 

(W(o=£^(o. 

a derivative of the approximation is *=° wherein ~Uk' 

represents the Fourier coefficients derivative . 

6. (Original) The method of Claim 5, wherein each coefficient u'k is 
computed from uk. 

7. (Previously presented) The method of Claim 5, wherein the circuit 
is a radio frequency (RF) circuit. 
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8. (Previously presented) The method of Claim 4, wherein the step of 
applying a collocation method comprises applying Chebyshev collocation to 
discretize the set of differential-algebraic equation. 

9. (Cancelled) 

10. (Previously presented) The method of Claim 19, wherein the set of 
differential-algebraic equations comprises at least one of: a set of initial-value 
differential-algebraic equations and a set of boundary-value differential-algebraic 
equations. 

1 1 . (Previously presented) The method of Claim 20, wherein the circuit 
simulation is a radio frequency (RF) circuit simulation. 

12. (Previously presented) The method of Claim 20, wherein the step 
of applying a collocation method comprises applying Chebyshev collocation to 
each differential-algebraic equation to discretize the set of differential-algebraic 
equations. 

13. (Cancelled) 

14. (Previously presented) The method of Claim 19, further comprising 
enforcing continuity of the solution at the boundary of neighboring intervals. 

15. (Currently Amended) The method of Claim 20, wherein the set of 
differential-algebraic equations comprises a set of boundary-value differential- 
algebraic equations discretized into intervals , and wherein the boundary-value 
differential-algebraic e quat i ons equation intervals include a first and a last 
interval. 
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16. (Original) The method of Claim 15, further comprising enforcing a 
boundary condition at a boundary of the first and the last interval. 

17. (Previously presented) The method of Claim 20, further comprising: 
solving the set of differential-algebraic equations using a Newton- 

Raphson iterative method; and 

in each Newton-Raphson step of the Newton-Raphson iterative method, 
solving a linear Jacobian system using a linear iterative method. 

18. (Cancelled) 

19. (Previously presented) A method of solving a set of differential- 
algebraic equations arising in a circuit simulation, the method comprising:, 

applying a collocation method to each differential-algebraic equation to 
discretize the set of differential-algebraic equations; 

forming a solution to the set of differential-algebraic equations based on 
the discretized differential-algebraic equation; and 

determining an order of accuracy desired in each interval; 

wherein: 

the set of differential-algebraic equations comprises a set of boundary- 
value differential-algebraic equations, and wherein the boundary-value 
differential-algebraic equations are discretized in intervals, and wherein 
neighboring intervals share a boundary; and 

the solution in a particular interval is smooth, and wherein the step of 
determining the order of accuracy desired in each interval comprises 
determining whether to increase the order of accuracy of the particular interval. 

20. (Previously presented) A method of solving a set of differential- 
algebraic equations arising in a circuit simulation, the method comprising: 
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applying a collocation method to each differential-algebraic equation to 
discretize the set of differential-algebraic equations; 

forming a solution to the set of differential-algebraic equations based on 
the discretized differential-algebraic equation; and 

determining an order of accuracy desired in each interval; 

wherein: 

the set of differential-algebraic equations comprises a set of boundary- 
value differential-algebraic equations, and wherein the boundary-value 
differential-algebraic equations are discretized in intervals, and wherein 
neighboring intervals share a boundary; and 

the solution in a particular interval is not smooth, and wherein the step of 
determining the order of accuracy desired in each interval comprises splitting 
the particular interval into at least two subintervals. 

21. (Original) The method of Claim 17, further comprising separately 
approximating for each interval a local preconditioner. 

22. (Original) The method of Claim 21, wherein the local preconditioner 
comprises at least one of: 

a capacitance matrix; and 
a conductance matrix. 

23. (Cancelled) 

24. (Previously presented) The computer-readable medium of Claim 34, 
wherein the set of differential-algebraic equations comprises at least one of: a 
set of initial-value differential-algebraic equations and a set of boundary-value 
differential-algebraic equations. 
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25. (Previously presented) The computer-readable medium of Claim 34, 
wherein the circuit simulation is a radio frequency (RF) circuit simulation. 

26. (Previously presented) The computer-readable medium of Claim 34, 
wherein the step of applying a collocation method further causes the processor 
to carry out the step applying Chebyshev collocation to each differential- 
algebraic equation to discretize the set of differential-algebraic equations. 

27. (Cancelled) 

28. (Previously presented) The computer-readable medium of Claim 34,, 
wherein the instructions further cause the processor to carry out the step of 
enforcing continuity of the solution at the boundary of neighboring intervals. 

29. (Currently Amended) The computer-readable medium of Claim 34, 
wherein the set of differential-algebraic equations comprises a set of boundary- 
value differential-algebraic equations discretized into intervals , and wherein the 
boundary-value differential-algebraic e quations equation intervals include a first 
and a last interval. 

30. (Original) The computer-readable medium of Claim 29, wherein the 
instructions further cause the processor to carry out the step of enforcing a 
boundary condition at a boundary of the first and the last interval. 

31. (Previously presented) The computer-readable medium of Claim 34, 
wherein the instructions further cause the processor to carry out the steps of: 

solving the set of differential-algebraic equations using a Newton- 
Raphson iterative method; and 

in each Newton-Raphson step of the Newton-Raphson iterative method, 
solving a linear Jacobian system using a linear iterative method. 
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32. (Cancelled) 

33. (Currently amended) A computer-readable medium carrying one or 
more sequences of one or more instructions for solving a set of differential- 
algebraic equations arising in a circuit simulation, the one or more sequences of 
one or more instructions including instructions which, when executed by one or 
more processors, cause the one or more processors to perform the steps of: 

applying a collocation method to each differential-algebraic equation to 
discretize the set of differential-algebraic equations; and 

forming a solution to the set of differential-algebraic equations based on 
the discretized differential-algebraic equation; 

wherein: 

the set of differential-algebraic equations comprises a set of boundary- 
value differential-algebraic equations, and wherein the boundary-value 
differential-algebraic equations are discretized in intervals, and wherein 
neighboring intervals share a boundary; 

the instructions further cause the processor to carry out the step of 
determining an order of accuracy desired in each interval; and 

the solution in a particular interval is smooth, and wherein the step of 
determining the order of accuracy desired in each interval further causes the 
processor to carry out the step of determining whether to increase order of 
accuracy of the particular interval. 

34. (Previously presented) A computer-readable medium carrying one 
or more sequences of one or more instructions for solving a set of differential- 
algebraic equations arising in a circuit simulation, the one or more sequences of 
one or more instructions including instructions which, when executed by one or 
more processors, cause the one or more processors to perform the steps of: 
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applying a collocation method to each differential-algebraic equation to 
discretize the set of differential-algebraic equations; and 

forming a solution to the set of differential-algebraic equations based on 
the discretized differential-algebraic equation; 

wherein: 

the set of differential-algebraic equations comprises a set of boundary- 
value differential-algebraic equations, and wherein the boundary-value 
differential-algebraic equations are discretized in intervals, and wherein 
neighboring intervals share a boundary; 

the instructions further cause the processor to carry out the step of 
determining an order of accuracy desired in each interval; and 

the solution in a particular interval is not smooth, and wherein the step of 
determining the order of accuracy desired in each interval further causes the 
processor to carry out the step of splitting the particular interval into at least 
two subintervals. 

35. (Original) The computer-readable medium of Claim 31, wherein the 
instructions further cause the processor to carry out the step of separately 
approximating for each interval a local preconditioner. 

36. (Original) The computer-readable medium of Claim 35, wherein the 
local preconditioner comprises at least one of: 

a capacitance matrix; and 
a conductance matrix. 

37. (Currently Amended) A method of simulating an rf circuit, 
comprising the steps of: 

determining a p l ura li ty plurality of differential equations describing 
operation of the rf circuit; 
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determining a set of Chebyshev Gauss-Lobatto collocation points for the 
plurality of differential equations, producing a set of intervals; 

discretizing each of the differential equations based on the Chebyshev 
Gauss-Lobatto collocation point intervals; 

determining a smoothness for each interval, increasing an order of a 
solution for an interval if it is smooth, and splitting the interval into at least two 
sub-intervals if the interval is not smooth; 

solving the differential equations in each of the intervals; and 

simulating the rf circuit based on the solved intervals. 

38. (Currently Amended) A The method according to Cla i m 37, of 
simulating an rf circuit, comprising the steps of: 

determining a plurality of differential equations describing operation of the 
rf circuit; 

determining a set of Chebyshev Gauss-Lobatto collocation points for the 
plurality of differential equations, producing a set of intervals; 

discretizing each of the differential equations based on the Chebyshev 
Gauss-Lobatto collocation point intervals; 

solving the differential equations in each of the intervals; and 

simulating the rf circuit based on the solved intervals; 

wherein the step of solving comprises applying a set of at least one high 
order solution to at least one of the intervals and applying at least one solution 
from a set of low order solutions to a plurality of the intervals. 

39. (Original) The method according to Claim [[37]] 38, 

wherein the step of solving comprises applying a set of solutions 
comprising more low order solutions than high order solutions. 

40. (Currently Amended) The method according to Claim [[37]] 38, 
further comprising the steps of: 
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dividing the intervals into smooth and non-smooth categories, 
applying higher order solutions to the smooth category intervals, and 
applying lower order solutions to the non-smooth category intervals. 

41. (Currently Amended) A The method according to C l aim 37, of 
simulating an rf circuit, comprising the steps of: 

determining a plurality of differential equations describing operation of the 
rf circuit; 

determining a set of Chebyshev Gauss-Lobatto collocation points for the 
plurality of differential equations, producing a set of intervals; 

discretizing each of the differential equations based on the Chebyshev 
Gauss-Lobatto collocation point intervals; 

solving the differential equations in each of the intervals; and 

simulating the rf circuit based on the solved intervals; 

wherein the Chebyshev Gauss-Lobatto collocation points produce a small 
number of intervals in areas in which the differential equations exhibit high 
convergence, and a large number of intervals in areas where the differential 
equations exhibit low convergence. . 

42. (Currently Amended) The method according to Claim [[37]] 38, 
wherein the step of solving comprises applying higher order solutions in smooth 
intervals and applying lower order solutions in less smooth intervals. 

43. (Previously presented) The method according to Claim 37, further 
comprising the steps of: 

enforcing continuity of the solution at each interval boundary; and 
enforcing a periodic boundary condition at each first and last interval 
boundaries. 
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44. (Previously presented) A method of simulating response in a circuit, 
comprising the steps of: 

determining a plurality of differential equations describing operation of the 
circuit; 

determining a set of collocation points for the plurality of differential 
equations, producing a set of intervals comprising at least one high 
convergence interval and a plurality of low convergence intervals; 

applying a higher order solution in the at least one high convergence 
interval; 

applying a lower order solution in the low convergence intervals; and 
simulating the circuit response using the higher and lower order solutions. 

45. (Previously presented) The method according to Claim 44, wherein 
the collocation points comprise Chebyshev Gauss-Lobatto collocation points. 

46. (Previously presented) The method according to Claim 44, further 
comprising the steps of: 

enforcing continuity of the solutions at each interval boundary; and 
enforcing a periodic boundary condition at each first and last interval 
boundaries. 

47. (Currently amended) The method according to Claim 44, further 
comprising the steps of: 

enforcing continuity of the solutions at each interval boundary; and 
enforcing a periodic boundary condition at at least one of the first and 
last interval boundaries. 

48. (Previously presented) The method according to Claim 44, wherein 
more low order solutions are applied to the differential equations than high 
order solutions. 
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